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Abstract 



We show that the statistical multifragmentation model with the standard param- 
eterization of the free volume predicts a constant pressure as density approaches the 
ON , normal nuclear density. This is contrary to the Raduta&Raduta result obtained by 

disregarding the center of mass constraint. It is demonstrated that in finite nuclear 
^r) ' systems the partitions with small number of fragments play an important role. 
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In recent years much efforts have been directed to investigating thermodynamical prop- 
erties of finite nuclear systems (see e. g. refs. [Ij]). In the paper |2| the authors study 
the phase diagram of finite systems with a statistical model designed for describing nuclear 
multifragmentation. The authors' model follows the theoretical prescriptions developed ear- 
lier in refs. [3-5] and based on the assumption that fragments are formed at a low density 
freeze-out stage. Consequently, the thermodynamical characteristics calculated within such 
kind of models have adequate physical meaning at low densities only, e.g. p<(| — |)po 
(po ~0.15/m -3 is the normal nuclear density), when individual fragments are singled out 
from the surrounding nuclear matter || ffl. However, the authors of ref. |2| have tried to 
apply their model for the higher densities p — > po by removing the physical assumption that 
fragments do not overlap. They have found that the model predicts a Van der Waals kind 
of phase diagram. In particular, the pressure P — > oo at p — ► p , as seen from Fig. 3 of 
ref. ||. In this comment we want to demonstrate that this result can be misleading, even 
under assumption of overlapping fragments. Within the standard statistical treatment of 
finite systems, when the conservation laws are properly implemented, the pressure always 
remains finite, even in the limit p — > po- 

In the statistical models [2-5] the volume (density) affects the partition probabilities via 
so called free volume Vf determining the phase space available for the translational motion of 
fragments. It differs from the actual physical volume of the system V because of the finite size 
of the fragments and fragment-fragment interaction. In the standard canonical description 

1 A comment on "Investigating the phase diagram of finite extensive and nonextensive systems" by Al.H. 
Raduta and Ad.R. Raduta, Phys. Rev. Lett. 87, 202701 (2001) 
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of multifragmental ion the free volume enters in the partition probabilities through the factor 
(Vf/X^) N ~ 1 , where N is the number of fragments in the partition and = 2irh/ ' y/m^T is 
the thermal wavelength of nucleons at temperature T. The exponent (N — 1) comes from 
the integration over momenta and coordinates of fragments under constraints that total 
momentum and center of mass position are fixed for all partitions. These constraints are 
crucial for finite systems in contrast to the thermodynamic limit (N — > oo , V — > oo , N/V = 
const), where this factor is usually taken as (Vf/\\) N . In the canonical ensemble, for a finite 
system with the number of nucleons Aq, the free energy can be represented as 

F — —T 




where Cat is the volume- independent weight of iV-fragment partitions which is influenced by 
combinatorial factors as well as by internal excitation of fragments and their interaction. For 
simplicity, below we disregard the Coulomb interaction in the system. In the microcanonical 
ensemble the structure of the statistical weights with respect to the volume will be practically 
the same, and thus our conclusions concerning the volume dependence of the pressure in the 
limit Vf — > will not change. Though the weights Cn can significantly vary from partition 
to partition, they are finite for any finite system. Therefore, the behavior of the free energy 
in the limit Vf — > is similar for different assumptions on the fragments' relative motion in 
the freeze-out volume. To illustrate our point we use the excluded volume approximation 
Vf — V — Vq, where Vq = A /p is the normal nuclear volume of all fragments. In this 
approximation one gets straightforwardly: 

dF _ EiUN - l)V f N ~ 2 C N 

ov EiUVf-'Cv • [Z) 

Now it is obvious that in the limit V — > V the pressure goes to a constant, P — > T- (C 2 /Ci) = 
const. It is interesting to note that a similar behavior for this kind of model was also found 
in the thermodynamic limit ||. In this case constant parts of pressure isotherms appear in 
the coexistence region of liquid and gaseous phases. And the liquid phase is represented by 
an infinite cluster which in a finite system would correspond to the compound-like nucleus. 

We believe that the main reason of the P — > oo behavior obtained in ref. |§ is in 
disregarding the center of mass constraint []. The unconstrained integration over coordinates 
of all fragments in the freeze-out volume results in a factor proportional to Vf . If in eqs. (0) 
and (H) we take Vj N instead of Vf -1 , we get formally P ~ 1/V/ at V — > Vo, i.e. P — > oo. One 
can easily see that Vf ree adopted in ref. (their eq. (4)) gives the same limiting behavior 
asV f N . 

For illustration, in Fig. 1 we show the phase diagram obtained with the canonical calcula- 
tions for partitioning a one-component system containing v4 =100 nucleons into all possible 
fragments characterized by the mass number A only. Contributions of all 190569292 par- 
titions were calculated directly, in oder to achieve the best accuracy [0. We adopt the 
liquid-drop description for individual fragments disregarding the Coulomb interaction and 

2 Because of this and other approximations used in the model M, depicting it as "exact" in comparison 
with other models does not seem justified. 
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Figure 1 : Pressure versus free volume (in units of normal nuclear volume) in multifragmentation 
of one component system with ylo=100 nucleons calculated within the canonical ensemble. Solid 
lines are calculations taking into account the center-of-mass constraint in the fragmentation sta- 
tistical weights, dashed lines correspond to a similar calculation but without this constraint. The 
temperatures (in MeV) are given by numbers at the curves. 



internal excitation of fragments. The statistical weight of a partition with N fragments with 
individual multiplicities Na is taken as 

N A 

(3) 

where B(A) = a v • A — a s • A 2 / 3 is the liquid-drop binding energy of fragment A with 
parameters ay ~ as ~ 16 MeV. In this case the pressure is P = T ■ ((N) — 1)/Vf, where (JV) 
is the mean fragment multiplicity. It is seen that at large Vf the phase diagram is consistent 
with expectations for a gas system. Here the system disintegrates into many fragments and 
the fragment mass distribution falls off nearly exponentially with A. By decreasing Vf we 
move into region where the partitions with low fragment multiplicity dominate. The mass 
distribution turns into "U-shape"-like one, which is associated with the phase transition. 



iV — 1 
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Namely in this region the pressure has a slight "backbending" (at low temperatures) and then 
approaches a constant. Formally this behavior at small volumes sets in because probabilities 
of the channels with N > 1 are suppressed by factors oc Vf N ~ l . The compound nucleus (a 
partition with N—l) dominates in this case, however, it does not contribute to the pressure. 
This is a trivial consequence of the conservation laws: the compound nucleus is at rest in 
its center of mass frame. For N > 1 only relative momenta and positions of fragments have 
physical meaning. 

We have also simulated the effect of non-conservation of the center-of-mass by introducing 
an additional factor oc Vf in the weights of all partitions in the formula (^]). This effect is 
shown by dashed lines in Fig. 1, corresponding to the pressure P = T ■ (N) /Vf. This means 
that compound nucleus can "move" and "exert a pressure" in the freeze-out volume, and 
as a result P — > oo at Vf — > 0. Moreover, in this case the phase diagram looks like a Van 
der Waals one, that in fact does not correspond to the physical content of the model. Even 
in the region of large volumes relevant for the model this effect can considerably change 
many thermodynamical characteristics, such as the critical temperature or the caloric curve 
at constant pressure. 

This example shows that phase diagrams of finite systems depend sensitively on physical 
assumptions adopted in the model. In particular, in the study of a liquid-gas type phase 
transition a careful treatment of the partitions with small multiplicities becomes extremely 
important. However, in order to investigate realistically the phase diagram at p > \p$ it 
is necessary to introduce a really new physics in the statistical models, e.g. instead of the 
picture of individual fragments surrounded by the nucleon gas to consider the bubbles of the 
nuclear gas inside the nuclear matter. 
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